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Abstract. We investigate the statistics of records in a random sequence {a’s(0) = 
0 ,Xb(1),--- ,Xb(ti) = Xb( 0) = 0} of n time steps. The sequence Xs(/c)’s represents 
the position at step k of a random walk ‘bridge’ of n steps that starts and ends at the 
origin. At each step, the increment of the position is a random jump drawn from a 
specified symmetric distribution. We study the statistics of records and record ages 
for such a bridge sequence, for different jump distributions. In absence of the bridge 
condition, i.e., for a free random walk sequence, the statistics of the number and ages of 
records exhibits a ‘strong’ universality for all n, i.e., they are completely independent 
of the jump distribution as long as the distribution is continuous. We show that the 
presence of the bridge constraint destroys this strong ‘all n’ universality. Nevertheless 
a ‘weaker’ universality still remains for large n, where we show that the record statistics 
depends on the jump distributions only through a single parameter 0 < /r < 2, known 
as the Levy index of the walk, but are insensitive to the other details of the jump 
distribution. We derive the most general results (for arbitrary jump distributions) 
wherever possible and also present two exactly solvable cases. We present numerical 
simulations that verify our analytical results. 
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During the last few years, there has been a growing interest in the study of records. 
Records have not only become popular in our societies as one often hears and reads, 
in the media, about record breaking events, including for instance sports [I] or weather 
records [2j, but they have also found interesting applications in various areas of 
sciences. Records are often useful to characterize the statistics of a (discrete) time 
series x(0),a;(l), ... ,x(n) where a record happens at time k if x(k) is larger than the 
previous values x(0), x(l), • • • ,x(k — 1). As such, they have been studied in several 
contexts ranging from domain wall dynamics [3], spin-glasses [1] and random walks 
0 El 0 El El ED! to avalanches HU, models of stock prices mmm, growing networks 
models [Qj or the study of global warming p~5i rroi R7] and also in evolutionary biology 
[18, 19] (see Ref. [20 for a recent review). 

When considering record statistics, the most natural questions are the following: 
(a) what is the number of records after time n? (b) what is the probability that a record 
is broken at step nl how long does a record survive and in particular what is the age of 
the longest lasting record? These questions are fully understood in the case where the 
variables x(*)’s are independent and identically distributed (i.i.d.) [2DI22J (see also Ref. 
[23 J for a short review), although more refined questions about first passage properties 
of records in this i.i.d. case were only recently investigated [211 G5J. However, many 
applications of records in statistical physics have highlighted the importance of strong 
correlations between the random variables x(i)’s, for which much less is known. With 
that perspective, it was demonstrated rather recently that one-dimensional random 
walks offer a very useful instance of strongly correlated variables where the impact of 
correlations on records can be studied analytically 0 El 0 El U ITU] . As we recall it 
below, a remarkable feature of the record statistics of RWs with continuous jumps is 
that it is completely universal, i.e., independent of the jump distributions, even for a 
finite number of steps [5j. It is thus natural to ask whether this universality still holds 
for constrained random walks. This question is the main motivation for the present 
work where we consider random walk bridges, where the walker is constrained to start 
and end at the same point, which are relevant, for instance, to study periodic correlated 
time series. 

Let us consider a time series where the x(k) , s are the positions of a random walker 
(RW) after step k, evolving according to the Markov rule: 


x(0) = 0 , x(k) = x{k — 1) + rj(k) , (1.1) 

where T]{k)' s are i.i.d. random variables. We denote by p{r]) their symmetric probability 
distribution function (PDF) [with p(j]) = p(—r))}. In the following, we will focus on 
the positions of RW bridges {xsik)}, 0 < k < n which are RWs as defined in (1.1) 
with the additional constraint that they return back to the origin after n time steps, 
x B {n) = x B (0) = 0 (see Fig. [If. 

For the purpose of our study, it is important to distinguish between jump 
distributions p{rj) that are respectively discrete and continuous. A representative of 
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P d (v) = \ S (V + !) + !) , (i- 2 ) 

where the superscript d refers to this particular discrete case. In this case, the walker 
moves on a lattice of the integers. In contrast, when the jump distribution is continuous, 
we will denote it generally by a superscript c, i.e., p c (r]). Let p c (k ) = dpp c (rj) e lkv 
denote the Fourier transform of the jump distribution. It turns out that for many 
properties associated with the record statistics, only the small fc-behavior of the jump 
distribution matters. Generically, for small k, one typically has 

p c (k) = 1 — \akf + o(|fc| M ) , (1.3) 


where 0 < p < 2 is called the Levy index of the walk and a > 0 just sets the length scale 
of the jump. For p — 2, the variance a 2 = (if) is finite. Consequently, a 2 = <t 2 /2 and 
the RW converges, for a large number of steps n, to the Brownian motion. In contrast, 
for p < 2, the jump distribution has a heavy tail, p c (rf) ~ for large r/, leading 

to a diverging second moment. A random walk with p < 2 is called a Levy flight with 
index p. Among the class of continuous jump distributions with a finite a 2 , a special 
role is played by the exponential distribution 

p e (v) = ( 1 - 4 ) 


where b > 0 sets the length scale of the jumps and the superscript ‘e’ refers to 
exponential. This clearly belongs to the class in Eq. (1.3) with p — 2 and a = b. 
We will see later that this exponential distribution belonging to the continuous family, 
and the lattice RW in Eq. (1.2) belonging to the discrete family, represent two rare 
cases where the record statistics for the bridge sequence is exactly solvable. 

When studying the sequence of records of a time series x(0), x(l), • • • , x(n), a record 
happens at step k if it is strictly larger than the previous entries, i.e., 


x(k) > max(x(0), x(l), • • • , x(k — 1)) 


(1.5) 


and by convention, a;(0) is a record. We emphasize that the strict inequality in Eq. 
( 1.5| ) is of special importance in the case of discrete random walks (1.2). We denote 
by R(n) > 1 the number of records after n steps. Other important observables for our 
study are the ages of these records (see Fig. [l]). We define t* as the number of steps 
between the fc-th and (k + l)-th records: this is the age of the fc-tli record, i.e. the 
time up to which the k- th record survives. Note that the last record occurring before 
n is still a record at step n and we denote by A n the current age of this record at step 
n (see Fig. [Tj) . The typical fluctuations of these sequences of ages is rather simple to 
study. For instance, the typical age of a record is simply related to the average number 
of records ( R(n)) via £ typ ~ n/(R(n)). However, it has been shown that the statistics 
of these sequences is dominated, for RWs, by rare events ©mg. To characterize the 
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x B {k) 



Figure 1 . A realization of a random walk bridge x B {k) with an exponential jump 
distribution (a = e) with n = 20 steps. Here the number of records is R B {n) = 5. The 
Ti s denote the ages of the records and the pd s are the increments between successive 
records. Finally, A n denotes the age of the last record. The joint PDF of the rfs, p,'s, 
A n and R e B {n) is given in Eq. (3.21) and constitutes one of the main result of this 
paper. 


fluctuations of these rare events, it is useful to consider the probability Q(n) that the 
last time interval A n is the longest one, which we call in the following the probability 
of record breaking. For a RW with a given number of records R{n ) = m, the sequence 
of the ages of the records is r 1; ..., r m _i and A n . It is then useful to introduce the joint 
probability Q(m,n ) defined as 

Q(m, n) = Pr(Ai > max(ri,..., r m _i), R(n) = m) . (1.6) 

The probability of record breaking Q(ri) is then obtained by summing Q(m,n ) over all 
possible number of records in the time series, i.e., 

OO 

Q(n) = . (1.7) 

m= 1 

Another important observable to characterize the rare and extreme fluctuations of the 
sequences of ages is the age of the longest lasting record ^ max (n) which is defined as 

4iax(n) = max(ri, r 2 ,..., r m _i, A n ) , (1.8) 


whose fluctuations, for RWs, were recently demonstrated to be very sensitive to the last 

record mu- 

Notations: In the following, we will denote by R a (n), Q a {n ) and f^axCn) the 


quantity defined above for free random walks (1.1) where the superscript a refers to 
different jump distributions (|1.2[ |1.3[ |1.4[), where 


{ d corresponding to discrete jump distribution p d (rj ) , 

c corresponding to generic continuous jump distribution p c {rf) , (1.9) 
e corresponding to exponential jump distribution p e {rf) , 
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and where p d (p), p c (r/) and p e (p) are defined respectively in Eq. (1.2), (1.3) and (1.4) 


Besides we will use the notations R%(n) 7 Q%(n) and 7'“ ax n( n ) f° r the number of records, 


the probability of record breaking (1.7) and the age of the longest lasting record (1.8) 


for random bridges, where the subscript B refers to bridges. 

Before turning to our results for such constrained RWs, let us first remind the main 
known results for the record statistics of free RWs, where the walker can end up at any 
position x n after n time steps Min]. We first focus on the case of continuous 

(symmetric) jump distributions p c (p ) (|1.3). In this case, a remarkable property is that 


the full statistics of records, including the number of records R c (n ) as well as Q c (n ) (1.7) 
and l ^ a . Y (n) (1.8) are completely universal, i.e., independent of the jump distribution 
p c fg) - including Levy flights - even for finite n [5]. This universal behavior stems from 
the universality of the Sparre Andersen theorem [26] • In particular, for large n, the 


( 1 . 10 ) 


average number of records ( R c {n )) behaves for large n as [5] 

^ 2 

( R c (n )) ~ A c y/n , A c = —= , as n —$■ oo , 

7 r 


independently of 0 < p < 2 characterizing the jump distribution p c (p) in (1.3). This 


behavior (1.10) should be compared to the logarithmic growth found for i.i.d. random 
variables 


Moreover, it is also possible to compute exactly the full probability 
distribution P c (m,n ) = Pr (R c (n) = m,n ) of the number of records which, for large n 
and large m keeping A" = m/y/n fixed, takes the scaling form [5] 


P c (m,n) ~ 

In 


TYl \ 1 

X = — ) where <p c (A) = —e" 
n / \/7r 


~r ,X > 0 . (1.11) 


For free random walks, it is also possible to compute exactly Q c (n ) in Eq. ( |1.7[ ) and 
show that it converges, for large n, to a non-trivial universal constant given by [5j 

1 


lim Q c 

n—>oo 


n 


= Q C 


oo = 


dx-- 


= 0.626508..., (1.12) 


qj 1 + y/n xe* eii(y/x) 
a constant which also arises in the excursion theory of Brownian motion m m- 
Interestingly, in the case of free RWs, this probability Q c (n ) turns out to be related 
to the average value of (f'jj 


pc 

"max V 




(n)} in Eq. (1.8) via the relation [27] 
c (n + !)) = (CaxM) + Q c (ri) , 


implying in particular that 

(4axW) 


lim 

n—>■ oo 


n 


= AU = Q c 


.oo 


= 0.626508 ... 


while the full distribution of /7i was studied more recently in 


(1.13) 

(1.14) 
(see also [30]). 


For discrete free random walks, with a jump distribution given in Eq. (1.2), it was 


shown in |5j that the statistics of records is quantitatively different. In particular, in 
this case, the mean number of records ( R d (n )) also grows like y/n but with a different 
prefactor 


(1.15) 


{R A {n)) ~ A A y/n , A d = — , as n —* oo , 
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Free RW 


RW Bridge 

jump 

distribution 

discrete 

a = d 

continuous 

a = c 


discrete 

a = d 

exponential 

a = e 

continuous 

a = c 

(R a (n)) 


2 r 

~ -7=Vn 
a/tt 

(R%(n)) 

\fk . — 

Z2 

0F , — 
-—■Jn 

~ X B {g)Vn 

Q a (n) 

~ Q d (oo) = Q c (oo) 

= 0.626508 ... 

Q %( n ) 

~ Qb( oo) — Qb( oo) 

= 0.6543037.. 

? 

(CaxM) 

~ A d = A° 

''max ''max 

= Q d ’ c (oc) 

(Cax,s( n )) 

~ A d — A e 

''max, 77 ''max, 77 

= 0.6380640.. 

? 

n 

n 


Table 1 . Summary of the results for the record statistics of RW bridges, obtained 
in the present paper. For comparison, we have also presented, in the left part of the 
table, the results for free RW, obtained in Refs. [S|. Note that, for the free RW, the 
results for continuous distributions are completely universal, and thus also hold for 
an exponential distribution, which is not true for t he RW bri dges. The expression of 
the constants Q B e (oo) and A d [ ax B is given in Eqs. (4.9) and (4.30) respectively. The 
computation of Q c B {n ) and of (f^ ax B (n)) f° r a continuous distribution with arbitrary 
Levy index g [see Eq. (1.3)] remains an open question - hence the question marks (?) 
in the table. 


which is 1 /a/ 2 of the expression for the mean in the continuous case (1.10). Similarly, 
the distribution of R d (n), P d (m,n) = Pr(P d (n) = m,n ), takes the 
form, for large n and large m keeping X = m/^/n fixed 00: 


’ollowing scaling 


P d (m, n) ~ -^=g) d (x = -^L) , (X) = J-e 

\/n V \/n / V vr 


, X > 0 , (1.16) 


which is simply obtained from its continuous counterpart tp c (X) in Eq. (1.11) with the 
substitution A" —> \[2X. On the other hand, the results for the probability of record 
breaking Q d {n ) behaves, for large n, exactly as in the continuous case (1.12), i.e. 

r°° i 

lim Q d (n) = </ d (oo) = Q c ( oo) = / dx- - — —-— = 0.626508 ... (1.17) 


1 + a/7 ix e x erf (a/S) 


Moreover, the relation in Eq. (1.13) also holds for the corresponding discrete quantities, 


from which it follows, using Eq. (1.17), that 

(")> 


lim 

n—> oo 




d 

max V 


n 


= C = Q (oo) = 0.626508 .... 


(1.18) 


For random walk bridges, the situation is quite different. As expected, the statistics 


of records for discrete (1.2) and continuous (1.3) jump distributions are still different, as 
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they are for free random walks. But in this case, for continuous distributions, the 
statistics of records, for finite n, is not universal any more and depends on p c {j]). 
One expects however, and we can show it explicitly for the average number of records 
(. R B (n )), that, for large n, the various observables characterizing the record statistics 


depend only on the exponent p (1.3) and not on the further microscopic details of the 


jump distribution p c (r)). Another important feature of the record statistics of random 
walk bridges is that it is technically much more involved. Indeed, for free random 
walks, the computations require the full joint distribution of the ages of the records 
Ti, t 2 , ■ • • , r m _ i, A n but there is no need to keep track of the actual value of the record 
at a given time step. The knowledge of the actual value of the record at a given time 
step is however required for bridges, where the random walk returns back to the origin 
after n time steps. This is done here by considering the full joint distribution of the ages 
Ti s and the record increments pfs (which are the differences between two consecutive 
records), see Fig. [Tj The computation of this joint distribution is the main technical 


achievement of the present paper [see Eq. (3.21) below]. Finally, yet another noticeable 
difference between free RWs and RWs bridges, from the point of view of records, is 
that there is no simple relation between Q B (n) and (^max.s( n )) f° r bridges, while they 
are directly related for free random walks ( |1.13 ). Hence, one expects that Q%(n) and 
((“ a xB {n)) generically lead to two different universal constants, which we can compute 
explicitly in the case of the discrete a = d and exponential a = e distributions. 

Our main results can be summarized as follows (see also Table [l|. First, for 


discrete random walk (1.2) we obtain exact results for the full distribution of the records 
P B (m,n) = Pr(R B (n) = m,n). In particular, for large n (and n even as a discrete 
random walk bridge has necessarily an even number of steps), we show that 

(Ri 




A d 

B 


n 


A d = xfL 
B 2 3 / 2 


(1.19) 


Hence it also grows like y/n, as for free RW (1.15), but a with a different prefactor: 
A d /A d = 7t/4 < 1. On the other hand, for large n and large m, keeping X = m/y/n 
fixed, the distribution P B (m,n ) takes the scaling form 


\/n 


m 


X = — ) , tp B (X) =4Xe 
n 


—2X 2 


X > 0 , (1.20) 


which is different from its counterpart for free RW (1.16). 


For continuous jump distribution p c (r)), we obtain the large n behavior of ( R c B (n )) 


as 


(a%M) ~ A%(n)y/n , as n —> oo , 


( 1 . 21 ) 


where A c B {p) is given in terms of an integral which depends explicitly on /i [see Eqs. (2.23) 


and (2.24)]. In the case n — 2, this amplitude can be evaluated explicitly as 


7r 


A b (p — 2) — ^ i 


( 1 . 22 ) 


which should be compared to the value of the amplitude A c = 2/y / 7r, independently 


of /i, for the free RWs (1.10), and hence A b (/i = 2 )/A c = 7t/4 < 1. By comparing 
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this result ( p2| ) with the one obtained for discrete random walks, we obtain that the 
ratio A d B /A c B (p = 2) = 1/^2, as for free RWs. Although the analysis of the statistics 
of R B (n), beyond the first moment, is quite difficult, it is possible to compute exactly 
the full distribution P B {m,n ) of the number of records R B (n ), for any finite n. for the 
case of the exponential distribution p e {jf) (1.4), which is representative of the case fi — 2 
In particular, for large n and large m, keeping X = mj\fn fixed, the 


[see Eq. (1.3) 
distribution P B (m,n ) takes the scaling form 


P B (m,n) - —<p% 

In 


m 


X = — \ , <p° B (X)=2Xe 


-x 2 


X > 0 


(1.23) 


which is different from its counterpart for free RW (1.11), and also slightly different 
from the corresponding distribution ip B (X) for the discrete bridges (1.20). 

On the other hand, for the record breaking probability Q%{n ), we obtain exact 
results for the discrete (a = d) and exponential (a = e) jump distributions. In 
particular, for large n, we show that these quantities converge to the same constant 
which can be computed exactly in terms of a single integral given in Eq. (4.9). This 
integral can be easily evaluated numerically with arbitrary accuracy, yielding 

lim Q B (n) = lim Q B (n ) = Q B (oo) = Q B ( oo) = 0.6543037... , (1-24) 


which is different from the one characterizing free RWs (1.12), Q B ( oo) > Q d (oo). 
Furthermore, for the discrete and exponential distributions, we compute exactly the 
average age of the longest lasting record R (n)) and show that 


lim 

n—>oo 




n 


= lim 

n—>oo 


(Cax^fa)) 


n 


= A°„ |B = AL, b = 0.6380640. 


(1.25) 


where, again, this constant can be expressed in terms of an integral given in Eq. (4.41). 
Note in particular that A" ax B ^ Q B ( oo), while the two corresponding quantities coincide 
for free RWs [see Eqs. (1.14), (1.18)]. Although we can not prove it, we expect that our 
results obtained for the exponential case for the large n behavior of the distribution of 
the number of records (1.23), the probability of the record breaking (1.24) and the age 
of the longest lasting record (1.25) are valid for any continuous distribution p c {jf) with 
finite a, i.e., with /i = 2. This conjecture is corroborated by our numerical simulations 
of the record statistics of RW bridges with a Gaussian jump distribution. 

The paper is organized as follows. In section [2j we compute exactly the 
mean number of records for RW bridges, ( R%{n )), for the discrete and exponential 
distributions (for any value of n in these cases) and for generic continuous distributions 
(in the large n limit). In section [3] we compute exactly the full distribution of the 
number of records R B (n ) for the discrete and exponential distributions for arbitrary n. 
Section [5] is devoted to the statistics of the ages of the records, where we compute Q%(n) 
and a g a i n f° r th e discrete and exponential jump distributions. In section [5] 

we compare our exact analytical results to numerical simulations before we conclude in 
section [6} Some technical details have been left in the appendices A, B and C. 
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Figure 2. A discrete random walk bridge of n = 20 steps. Here the number of records 

is i?20 = •Emax,i?(20) + 1=6. 


2. Mean number of records 


Let R%(n) be the number of records up to step n for a random walk bridge. We can 
then write 


R“ B ( n ) = ^° a (k,n) , 


( 2 . 1 ) 


k =0 


where a a {k,n) is a binary random variable taking values 0 or 1. The variable 
cr a (k 1 n ) = 1 if a record happens at step k and a a (k,n ) = 0 otherwise, for a random 
walker which comes back to the origin after n steps. By convention a a (k = 0, n) = 1 
as +_b( 0) = 0 is a record. On the other hand one has obviously cr a (k,n ) = 0 for k > n 
[in particular a Q (?r,n) = 0 because xg(n) = xb( 0) = 0 is never a record - as a record is 
defined by a strict inequality in Eq. ( |1.5 )]. Hence, the mean number of records up to 
step n is given by: 


71—1 


71—1 


(R%(n)) = a (k,n)) = 5T“+ n ) , 


( 2 , 2 ) 

k =0 k =0 

where r a (k, n ) = (cr a (k, n )} is the record rate, i.e., the probability that a record happens 
at step k. Therefore, to compute (R^(n)), we will Erst evaluate r a (k,n) and then sum 
over k (2.2). For n > 1, the evaluation of r a (k,n ) relies on the two following quantities: 


• The free Green’s function (propagator) G a (x,xo,n) that denotes the probability 
(discrete distribution) or probability density (continuous distribution) that a 
random walker starting at +o arrives at x after n steps. 

• The constrained Green’s function Gf (x, xo, n ) that denotes the probability (discrete 
distribution) or probability density (continuous distribution) that a random walker 
starting at x 0 arrives at x after n steps and staying strictly positive in-between. 


To compute r a (k,n ), let us suppose that a record happens at step k with a record 
value x (see Fig. [2]). This corresponds to the event that the walker, starting at the 
origin at step 0, has reached the level x for the first time at step k and returns back 
to the origin after n steps - as we are considering random walk bridges. In the time 
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interval [0, k\, the walker propagates from 0 to x, being constrained to stay strictly below 
x. To compute the corresponding propagator, we take x as the new origin of space and 
then reverse both the time and coordinate axis. Hence, we see that on the time interval 
[0, k], the particle propagates with G" (x, 0, k ). On the other hand, between step k and 
step n (where the walker ends at the origin) the walker is free and thus propagates with 
G"(0, x, n — k ) = G"(x, 0, n — k), as the jump distribution is symmetric. The record rate 
is then obtained by integrating the probability of this event over x > 0 as the record 
can take place at any level x > 0 (note that only the first record, i.e., k = 0, is such 
that x = 0). Using the statistical independence of the random walk in the time intervals 
[0, k) and [k,n] (being Markovian), one thus has, for n > 1: 


r a (k,n ) = 


G“(0, 0, n) 


dxG"(x, 0,fc)G Q (x, 0,n — k) , 0 < k < n — 1 , (2.3) 


where we have divided by G“(0,0,n) as we are considering random walks that are 
conditioned to come back to the origin after n time steps (bridges). Note that in the 
case of a discrete random walk, i.e., a = d, the integral over x in Eq. (2.3) has to be 
replaced by a discrete sum. Let us now analyze this formula (2.3) for different types of 
random walks. 


2.1. Discrete random walk 


In the case of a discrete random walk, it is convenient to compute (-Rg(n)) = 
ELo rd (M) as 

(RbH) { o) 




(2.4) 


G d (0, 0, n) 

where, here and in the following, the subscript ‘(0)’ refers to a random walk starting 
from x(0) = 0 and ending at x(n) = 0 after n steps and G d (0,0,n) is the probability 
for this event. The generating function of the numerator (i?^(n))(o) is from Eq. (2.3) 


^V(/4(n))( 0 ) = ^G d (x,0,z)G d (x,0,z) , 


n =0 


x=0 


where the propagators can be explicitly computed (see Appendix A): 

1 


G d (x, 0, z) = f>"G d (x,0 ,n) = ^== ( - \ 

71=0 V \ / 


and 


a d > (x,0,z) = J^z"GUx,0, 


71 = 


71=0 


l-vr 


(2.5) 


( 2 , 6 ) 


(2.7) 


Hence one obtains by summing over x in Eq. (2.5), using Eqs. (2.6) and (2.7) 

1 


J2 z "( Rd B( n »( 


o) 


+ 


71=0 


2(1 ~z 2 ) 


( 2 . 8 ) 
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Therefore one obtains 


<tf fl(n )) ( „, ={ \ + 2^1® ’ f ° r " eV6n ’ n = 2k 

0 , for n odd 


(2.9) 


Using that the denominator in Eq. (2.4) is given by G d (0,0,n) = 2 n ( 1 f 2 ) if n is even 
and 0 otherwise one gets finally 


0 , n odd . 


( 2 . 10 ) 


The first terms are {R%( 0)) = 1, {R%( 2)) = 3/2, (R%( 4)) = 11/6. Finally, for large n it 
behaves like 

<flj(n))~{ tOIn eVen ' (2.11) 


0 , for n odd , 


as announced in Eq. (1.19) 


2.2. Exponential jump distribution 

In this case, the jumps r/j’s are continuous variables with a distribution given by 
P e ( 7 l) = As we will see, the record statistics for the exponential distribution 

is exactly solvable. Indeed, in this case, both the free and the constrained Green’s 
functions can be computed exactly [131] • The associated generating functions (GFs) 
G e (x, 0, z) (free) and G>(a;, 0,z) (constrained) read [31J (see also Appendix B) 


OO 


G e (x, 0, z) = > z n G e (x, 0,n) = 


^ 2 by/l 

n=1 
oo 


= e b 


- ^Ty’ClIi'.O. n) = 6(x) + 1 


( 2 . 12 ) 


(2.13) 


n =0 


Note that the definition of the GF of the free propagator G e (x,0,z) in Eq. (2.12) does 


not contain the n — 0 term - as it will not enter into the calculations - while this term 
n = 0 has to be retained for the constrained propagator and produces a delta function. 
The average number of records is given by ( R%(n )) = Ylk=o re (k,n) where the 


records rates r e (k,n ) are given by Eq. (2.3). To proceed, we write as above in Eq. (2.4): 


G « (0j0jn) ,n>l, 

where the numerator (i?g(n))( 0 ) is given by 

n ~ 1 poo 

(R e B (n))( o) = / dx G>(x, 0, k)G e (x, 0, n — k) , n > 1. 


k =0 


Hence its GF is given by 

OO 

^ 2 "<fl|(n)) (0) = 


dx G>(x, 0, z)G e (x, 0, z) = — 


n= 1 


+ 


4 b V1 - 2 v 1 — ;; 


(2.14) 


(2.15) 


,(2-16) 
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where in the last equality we have used the explicit expressions of G e (x, 0,z) and 


G>(x, 0, z) given in Eqs. (2.12) and (2.13). From Eq. (2.16) we easily extract (R e B (n))( 0 ) 


for all n. Finally, from Eq. (2.14) one obtains the mean number of records as 


i cy2n—3 

(R e B m = 1 , (-Rfl(w)) = ~ + 721P2Y for 71 > 1 > 

V n —1 ) 


(2.17) 


where we have used G e (0, 0, n) = ( 2n _ 2 )/(62 2ri x ), for n > 1. Interestingly, by comparing 


Eqs. (]2.10p and ( |2.17[ ) one sees that ( R e B (n )) = (R B (2n — 2)). For large n, one hnds 

(2.18) 


straightforwardly from Eq. (2.17) 


(R B (n)} ~ — Vn , 


which is \/2 times the result obtained for the discrete random walk in Eq. (2.11). 


2.3. Stable jump distribution 

Now we consider the more general situation where the jumps are continuous random 
variables distributed according to p c (r /) as given in Eq. (1.3). Although an exact 


calculation of ( R B (n )) for any finite n seems quite difficult in this case, one can perform 
a large n asymptotic analysis as follows. 

We recall that the average number of records is given by (. R c B (n )) = Ylk=o rC (^^ n )- 
One can show that this sum over k is dominated by the values of k ~ 0(n) which are 


thus large, when n 1. Hence, to evaluate the record rates r c (k,n) given in Eq. (2.3) 


for large k one can replace the propagators G c (x , 0, n—k ) and G>(x, 0, k ) by their scaling 
forms valid for k,n 1, with k/n fixed, and x 1, with x/n 1 ^ fixed. One has indeed 

1 


G c (x , 0 ,n — k) 
Gf(x, 0, k) 


-R 


x 


a^n — k) 1 /^ \a(n — kyR 

1 -r' x 


a^/nk 1 / 2 -^ 1 /v + Kak 1 /^ 

where the scaling functions are normalized, i.e., dx R(x) = 1 and J 0 °° dx R^ 
The scaling function R(x) is a Levy stable distribution: 

R( x ) = — dke~ ikx e~ w \ 

2?T I _ 


[x 


(2.19) 

( 2 . 20 ) 
_ = 1. 

( 2 . 21 ) 


and in particular -R(O) = T(1 + 1 //x)/7r. On the other hand, there is no explicit 
expression for R+(x) for generic fi < 2, while for p = 2 one has R+(x) = 2xe~ x2 . 

With such a normalization one can check in particular that by integrating 


G>(x,0,k) in Eq. (2.20) over x one recovers the survival probability q(k), which is 


the probability that the walker, starting from x(0) = 0 stays positive up to step k: 

1 


dx G>(x, 0, k ) = q(k) 


Jo y/n k 

in agreement with the Sparre Andersen theorem 


, as k —>■ oo , 


( 2 . 22 ) 


By inserting these scaling forms 
(2.19, 2.20) into the expression for r c (k,n ) in (2.3) one folds, that for large k and n 
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keeping k/n = y, fixed (with 0 < y < 1): 

r c (k,n) = —i=H (y — - 


n 


n 


H(y) = 


7T 


r(i + i//i) ^y(i- y y/^ J 0 


dx R + (x)R 


x 


(: y~ l - i) 1/V 


Finally, from this scaling form for the record rate (2.23), one obtains finally 

fi 


(2.23) 


{R c B {n)) = r c (k,n ) ~ A c B (y)^Jn , A c B (y) = / d yH(y) . (2.24) 

k =0 


Hence ( R B (n )) also grows like \fn for bridges but with an amplitude which depends on //, 
as announced in Eq. (1.21). In particular, one can check that A c B (y = 2) = i/T/2, which 
coincides, as expected, with the result obtained in the exponential case, see Eq. (2.18). 


3. Full distribution of the number of records 

In the previous section, we computed the average number of records ( R%{n )). Here we 
show that, in some cases, namely for discrete (a = d) and exponential (a = e) jump 
distributions, the full distribution of R%(n) can be computed exactly, for any finite n. 


3.1. Discrete random walks 


In the case of discrete random walks, the distribution of R%{n) can be computed by 
using the close relation between R B (n) and the maximum of the random walk up to 
step n, denoted by x maXi B(n) = max 0 < m <„XB(m) [8] (see also Fig. [2]): 

R%(n) = x maXjS (n) + 1 . (3.1) 


To derive this relation (3.1), let us consider the time evolution of the two processes 
Rs(n) and x maXi £(?7,). At the next time step n + 1, if a new site on the positive axis 
is visited for the first time, the process x maX) s(?r) increases by 1, otherwise its value 
remains the unchanged. On the other hand, when this event happens, then the record 
number R B (n) is also increased by 1, and otherwise it remains unchanged. Therefore we 
see that the two processes are locked with each other at all steps: for any realization of 
the random walk, we thus have x maXj £(r7 + 1) — x maX) B(?7,) = R%{n + 1) — R B (n). Given 
that, initially one has X m ax r(0) = xb(0) = 0 and i?^(0) = 1 (since, by convention, the 
first position is a record) one immediately obtains the relation in Eq. (3.1). Hence, this 
relation allows us to compute the PDF of R B (n) as 


P B {m,n) = Pr (R%{n) = m) = Pr(x maXjB (n) = m - 1) (3.2) 

= ^ , Pj(m,n)(o) = - !,0 ,fc)G|(m — 1,0 ,n - k) (3.3) 

where G>(x,x 0 , k) is the probability that the random walker, starting at x 0 , arrives at 
x after k steps, while staying non-negative (i.e., it may touch 0 but not —1) in-between. 
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^F^(m,n) (0 )*" = G>(m - 1,0,z)G>(?n - 1,0, z) , 


(3.4) 


n =0 


where G d (m — 1, 0, z ) is given in Eq. (2.7) and G d (;r, 0, z ) is given by (see Appendix A) 

x+l 


Gi{x,0,z) = ^Gi{x,0,n)z n = ( 


2 /I - vT^" 


n =0 


Therefore Eq. (3.4) reads explicitly 


J2 P B( m ^n)( 0) z n = - 


2(1 — \/T 


n=0 


2m—1 


(3.5) 


(3,6) 


Here, to perform this computation, we show an alternative method which will also 
be useful to compute the statistics of the ages. To do so, we first introduce the ages of 
the records r (s and write the joint probability distribution of the r (s (see Fig. [Tj) and 
the number of records R%(n), 

P d (£ l, • • • , 4i-i, a, m, n ) = Pr(ri = ..., r m _i = £ m _ i, A n = a, P|(n) = m, n) 

P d (^m,?i) ( o) 


G d (0,0,n) ’ 

where the numerator P d (£, m, n)(o) is given by 


^ m—1 


P d (^ m, n) (0) = / d (^i) ■ ■ • / d (C-i)G|(m - 1, 0, a)5 + a, 


n 


(3.7) 


(3,8) 


. i= 1 


where f d (£) is the first passage probability that the discrete RW, starting from xq, arrives 
at x 0 + 1 for the first time at step £, and S(a,b ) denotes a delta Kronecker function. 
Since the RW is invariant under translation, this probability is independent of xq and 
for a discrete RW, its generating function is given by 

1 - 


/ d M = £/■*(<) 


/ = 


(3.9) 


t=\ 

from which we deduce that 


f\£) = 


0 , £ even , 
(_!)(■«— 1)/2 


7T 


2T(1 — f/2)T(3/2 + £/2) 


,£ odd . 


(3.10) 


This joint distribution (3.7) contains all the information about the observables that 
we wish to compute here. In particular, the probability distribution of the record number 
R%(n) is obtained by integrating it over the £(s: 


OO OO 


P|(m,n) = ^ ^ ... Y pd ^ £ ^'" Jm-i,a,rn,n) = 

4 = 14=1 e m -i=ia,=i 


P^(m, 77.) (Q) 

G d (0, 0, n) 


, (3.11) 
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where the generating function of the numerator is given by 

OO 

^Vp|(m,n) ( o) = [f d (z)] m - l Gi{m - 1,0, z) 


n =0 


2 (1 - vT^" 


2m— 1 


(3.12) 


where we have used Eqs. (3.5) and (3.9). Thus we see that this result (3.12) coincides 


with the one obtained previously by a quite different method (3.6). From this expression 
(3.12) we can compute P|j(m, n)( 0 ) (using Cauchy’s formula) and finally Pg(m, n) as (for 


n even 


P a (m n) = 2" +1 (—1)-/ 2 +- 2 yt /2m - 


F(j/2 + 1) 


T(j/2 + 1 — n/ 2 — m) 


, (3.13) 


for 1 < m < n/ 2 + 1 (while P^(m,n) = 0 for m > n/ 2 + 1, as the maximum value 
of a discrete RW bridge can not exceed nj 2). Although the moments (Rg(n)) can be 


obtained, in principle, from the full distribution in Eq. (3.13) it is more convenient 


to compute them from the GF in (3.12). In particular, one can easily check that one 
recovers the result for the average number of records (P^(n)) as obtained previously 


(2.4), (2.8). Besides one can compute higher moments of Rg(n). For instance, the 
second moment is given by 


(Bw]A 

For large n, it behaves as 

<[/?!(«)] 2 >~ 


n+ 1 yV r(n/2 + 1) 

2 2~ T(n/2 + 1/2) ’ n ° dd ’ 

0 , n even . 


^ , n odd , 
0 , n even . 


(3.14) 


(3.15) 


Finally, from Eq. (3.6), one can also extract the expression of (rn, n) in the large n 


limit, which can be obtained by studying the limit z —> 1 of the GF in Eq. (3.6). In this 
limit, we set z = e -s and investigate the limit s —>• 0 where the discrete sum defining 
the GF in Eq. (3.6) can be replaced by an integral. Recalling that Pf,(m,n) = 0 if n is 


odd, the left hand side of Eq. (3.6) can be written, in the limit s —» 0 


J2 P B(m,n ) {0) e sn = ^ p b ( 2 k , ™) ( o) e 


— 2 sk 


n =0 


k =0 

1 / 

2 /,. 


d»e , ‘PB( m <y)(o ) ■ 


(3.16) 


can be written 


On the other hand, the right hand side of Eq. (3.6) assumes a simpler form in the limit 
z = e -s —> 1 such that Eq. 

I 


dye 


-<“PHm,y) ( o ) = 2e- 2 ”' /51 . 


0 


(3.17) 
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Hence Pg(m,i/)(o) can be straightforwardly obtained by Laplace inversion as 


Pg(m,n)( 0 ) 


J 2 - 


m 


7T n 


3/2 


r 2m2/n , for n even . 


(3.18) 


Finally, using that G d (0, 0, n) ~ ^/2/7Tn 1//2 , one hnds that the distribution P B {m,n) 
takes the scaling form (valid for n even): 


n 


P d (m,n) ~ 4=^ d ( X = ) , <p%(X) = 4Xe 


-2A' J 


(3.19) 


as announced in the introduction in Eq. (1.20). Note that this distribution ip B (X) 


coincides, up to a scale factor, with the probability distribution function of the maximum 
of a Brownian bridge on the unit time interval, as expected from the relation stated in 


Eq. (3.1). In particular, at variance with the result for the free random walk (1.16), 


this PDF is not a half-Gaussian distribution. In addition, it is easy to check that, from 
the moment of this distribution (p B (X), one re cover s the asymp totic result for (R B (n)) 
and ( [R b (n)] ) obtained respectively in Eqs. (12. lit) and (13.15h. 


3.2. Exponential distribution 


For the exponential jump distribution p e {rj) (1.4), the starting point of our analysis is 


the equivalent of the joint distribution given, for discrete random walks, in Eq. (3.7). 


However, because p e (r /) is a continuous distribution, this computation is more delicate 
than in the discrete case. Indeed, as we are considering random walk bridges, the weight 
of the last part of the paths, where the walker comes back to origin, i.e., the last segment 
of duration A n = a (see Fig. |lT) , involves the propagator G>(Y, 0,a) = G>(Y, 0, a) 


[see Eq. (3.3)] where Y = x maX) s(n) is the actual value of the last record, which 
coincides with the maximum. For discrete random walk the number of records R%(n) 
and x maX; s(n) are directly related through i mm jR (n) = R B (n) — 1 but this relation does 
not hold for a continuous jump distribution. Consequently, we need to keep track both 
of the number of records and of the value of the last record. A convenient way to do 
this is to introduce the record increments pf s and the joint distribution of the rf s, pf s, 
R e B (n) and x max , B (n) (see Fig. |Tj) : 

Pr[{rj = 4 pi G 4, n + dri]}i<j< m -i, A n = a, R e B [n) = m, x m ^ B (n) e [Y, Y + d Y},n\ 
= p e m Y)dr 1 ■ ■ ■ dr m _idE , (3.20) 

where the joint PDF P e ({4 a, m, Y ) is given by 

^ 771—1 „rv-i 


P e {{Zi,ri}i<i<m-i,a,m : Y) = 


G e (0, 0, n) 


n 

i=1 


d yi GUyi, 0, U - 1 )p e {yi + r t 


( m—1 


( m—1 


x G^(y,0,a)5 f J2 r i* Y ) 6 ( 5^4 + a,n J ,(3.21) 

where we have used that, for the exponential jump distribution which is continuous, 
G%(x,0,n) = Gt > (x,0,n) for x > 0. In particular, the joint distribution of the rf s, 














Record statistics for random walk bridges 


17 


A n and R e B (n), i.e., the equivalent of Eq. (3.7) for the discrete case, is obtained by 
integrating the formula in ( |3.20| ) over rf s and Y: 

P e {£i,£ 2 , ■ ■ ■ ,£m-i,a,m,n) = Pr(ri = £ 1: ..., r m _i = £ m -i,A n = a, R%(n) = m,n ) 

_ P e (£,m,n) { o) 


G e (0,0, n) ’ 


(3.22) 


where 


m-i roo 


IIO -L n, 

P e (£, m, n) (0 ) = / 

i=i • /o 


dr i T(£i, rf) 


x 


5 A + a,n\ , 

(3.23) 

r) . 

(3.24) 


where 

/»oo 

T(£,r)= dyG>(y,0,£-l)p e (y + r) . 

■Jo 

Note that this formula ( J3.23 ) together with (3.24) is actually valid for any continuous 
jump distribution p c (rf) (1.3)- where the superscript ‘e’ is replaced by ‘c’. However, 
its analysis is in general very hard to do, mainly because the constrained propagator 
Cr>(x, 0, n) does not have any explicit expression, which prevents one to perform the 
analysis of this multiple integral. Fortunately, such an explicit expression exists for the 
case of an exponential jump distribution p e fg) = 1/(2 b)e~^ v ^ b , which we now focus on. 

In this case, the generating function of the constrained propagator G>(x,0,?r) is 
given by Eq. (2.13), from which one gets that the GF of the building block T(£,r) in 
Eq. Q3.24[ ) is given by [using Eq. ( |2.13D 




z' — z 


i =i 


dyG e > (y,0,z)p e (y + r) = -(1 - y/l - z)e r/b . (3.25) 


From Eq. (3.25) one obtains 


T(l,r) = T ni)e- r,b , J2f‘ 


= i - Vi 


1=1 


which yields the expression of the coefficients / e (£) as 

At 1 


m = (-i) 


e+i 


, as t —> oo . 


(3.26) 


(3.27) 


m 


2 T(3/2 — £)Y{£ + 1) 2^ 3 / 2 

By comparing this expression for f e (£) and the one for f d (£) in Eq. (3.27) we easily see 
that / e (l) = / d ( 2£ — 1). By injecting this explicit expression of T(£,r) (3.26, 3.27) 
Eq. (3.23), the joint probability distribution P e (£, m,n)( 0 ) can be written 

m — 1 /.qo 

P e (l m, n) ( o) = J] f e {£i) / dYGl(Y, 0, a) e~ Y / b 


i =1 


X 


m_1 r°° Ay ( rn X \ f m ~ l \ 

n/ -j 8 w^ ri ~ y r (E^+ a ’ n ) 

i =1 J0 \ 2=1 /\ i =1 / 


(3.28) 
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Finally, using the identity 

m— 1 


m -1 noo ( m ~ 1 \ 

I] l dr^f^ri-yj 


Y 


m— 2 


(3.29) 


(m — 2)! 

which can be easily shown by taking the Laplace transform on both sides of ( 3.29[ ) with 
respect to Y, we obtain an expression for the joint probability of the rf s, A n and R%{n) 
as 

771—1 / m —1 

P e (f, m, n)( 0 ) = ]^[ f e (£i)q e (m, a) 5 ( ^^i + a,n ) , (3.30) 


2=1 

1 


. 2=1 


g (ra, aj = 


(m - 2)\b m ~ 1 J 0 


d y e - Y/6 y m - 2 G^(y,o,a), 


(3.31) 


which has thus a structure very similar to the one found in the discrete case (3.8), but 


with different building blocks. Furthermore, the GF of q e (m,a ) in Eq. (3.31) can be 
obtained explicitly as (see |Appendix B) 


q e (m, z) = ^ g e (m, a) z a = 
< 2=1 


1 1-yT^ _ (1 - y r Y^~zy 

b (1 + a/1 — z) m ~ l bz 


m—1 


. (3.32) 


From this joint distribution (3.30), together with Eqs. (3.26) and (3.32), it is 


possible to compute the statistics of all the observables which we want to study in this 
paper. In particular, one obtains the distribution of the number of records P^(m, n) by 


summing it over G,^ 2 , • • • Am -1 and a as done before in the discrete case (3.11). This 
yields 


P%{rn,n) = 


P|(m,n) ( o) 
G e (0, 0, n) 


^z n P e B (m,n) { o) = [f e {z)} m 1 q c (m, z) = J 

n =0 


1 (1-vT^i)' 


b (1 + Vl - ^) m_1 
1 (1 - ^/Y^~zf m - x 


y m— 1 


from which one gets 

o2n—1 f _1 2m—1 

P%(m,n) = m V (-1V 

( ) (n + m- 1 ! ^ V ; 

V n—1) v ' j =o 


2 m — 1\ 


r(j /2 + i) 


j J T(2 + j/2 — n — m) 


(3.33) 


(3.34) 


for 1 < m < n, which is independent of the scale parameter b [the factor 1/6 in 


Eq. (3.33) is indeed cancelled by the denominator in the first line of Eq. (3.33), 


G e (0,0, n) = (^Zf)/ (b 2 2n 1 )]. By comparing this result for the exponential distribution 


(3.34) with the corresponding one for discrete jumps obtained before in Eq. (3.13), we 
easily find that P e (m,n ) = P d (m,2n — 2), for n > 1. The moments of R%{n ) can be 


obtained, in principle, from the explicit expression of the full distribution in (3.34). It 


is however simpler to compute them from the GF in Eq. (3.33). For instance, from this 


expression, one can easily recover the value for ( R%(n )) obtained by the previous method 
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given in Eq. (2.17). Furthermore, from this relation (3.33), one can also compute higher 


moments of R B (n). In particular, the second moment is given by 


X>"([flBM] 2 }(0) 


+ 


n=0 


46(1 — z) 3 / 2 46(1 — z) 


from which one obtains, for n > 1 

l2 


([i^(n)] z ) =n--+ \^_[>) 
u BK n 7 2 2 T(n- 1/2) 


~ n , as n —> oo . 


(3.35) 


(3.36) 


Finally, from Eq. (|3.33|), one can obtain the limiting scaling form of the distribution of 

fixed, as 

, X > 0 , 


the number of records P B (m,n) for m 1, 1, keeping nnj\fn fixed, as 

m \ a „ _ X 2 


P%(m, n) ~ -W' fl [ x = ~7= ) , vi(X) = 2X e' 


n 


(3.37) 


where we have used that G e (0,0,n) ~ l/(26y / 7m), for n>l. This formula (3.37) yields 
the result announced in the introduction in Eq. (|1.23 ). Here also, one can easily check 


that, from the moments of the distribution tp e B (X) in Eq. (3.37), one recovers the large 
n behavior of (. R e B (n)) and ([i?£(n)] 2 ) obtained respectively in Eqs. (2.18) and (3.36). 


4. Statistics of the ages of records 

In this section, we study in detail two quantities associated to the ages of the records: 


(i) the probability Q%(n) (1.7) that the age of the last record is the longest one for a 
random walk bridge and (ii) the statistics of ax R (n) (1.8), which is the age of the 


longest lasting record. As in the previous section, we focus on the discrete (a = d) and 
the exponential (a = e) distributions. 


4.1. Probability of record breaking Q%{n) 

The computation of Q%(n) for a random walk bridge and for arbitrary jump distribution 
p(rf) is in general a very hard task. We show here that it can be computed exactly for 
the two special cases studied above: the discrete random walk (a = d) and the case of 
exponential jump distribution (a = e). 


4-1.1. Discrete random walk. The probability Q B (n ) 
from the full joint PDF of the intervals given in Eq. 


(1.7) can be simply computed 

(3.7) by summing it over the 
appropriate variables. For given values of the number of records m and of A n = a, we 
have to sum the joint PDF over the variables ifs from 1 to a - as A n = a is the longest 
interval (see Fig. [l|. Finally, we sum over all the possible values of a > 1 over the 
number of records, for m > 1. 


Q%(n) = pd (l m i n ) > 

m =1 a= 1 q = l £ m _ 1= l 


(4.1) 
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where P d (£,m,n ) is explicitly given in Eq. (3.7) and (3.8). Again, we separate the 
numerator and the denominator and write 

0 d , n) = QiWm 

G d (0, 0,n) 

oo oo a 

E^bWm*" = E : “E [h d (z,a)] x Gi(x,0,a) , 

n= 0 a= 1 cc=0 

where we have made the change of variable x = m — 1 and where 


(4.2) 


h d (z,a) = EZ / d ( fc )' 


(4.3) 


k =1 


From Eq. (4.2) one can obtain, in principle, the value of Q%{n)(Q) using Eq. (4.3) 
together with the explicit expression of (x, 0, a) given in Eq. (A.5) - for instance using 
Mathematica - though obtaining a closed form expression for Qg(n)( o) for any n seems 
quite difficult. It is however possible to extract the large n asymptotic behavior Q%(n)( 0 ) 
by analyzing Eq. (4.2) in the limit z = e -s when z —> 1, which thus corresponds to 
s —y 0. In this limit, the double sum on the right hand side of Eq. (4.2) is dominated by 
large a and large x. In this limit, keeping x/y/a fixed, G>(;r, 0,a) admits the following 
scaling form (for x + a even): 



2 1 


G > (x, 0, a) ~ 2\j -—-j=g \ — ),g(y) = ye 


X 


y 

2 


while G>(:r,0,a) = 0 for x + a odd. Similarly, using the expression of f c 
Eq. (3.10), we obtain that in the scaling limit a —> oo, s —> 0 (i.e., z = e“ 
the product sa fixed, h d (z,a ) in Eq. (4.3) takes the scaling form 

r d u .. 


h d (z,a ) ~ 1 — V2sF(as) , F(y) = 1 + 


2 A J 

= erf (y/y) + 


y 

1 e~y 

Vv 


(4.4) 

given in 
—> 1), keeping 

(4.5) 


From these asymptotic behaviors in Eqs. (4.4) and ( |4.5[ ) one obtains that the sum over 
x in Eq. (4.2), which becomes an integral in the limit s —> 0, a —» oo keeping y = as 
fixed, takes the scaling form (as a function of the scaling variable y) 

Td («> 


^2 [h d (z,a)] x G>(x,0,a) ~ 2 - V2sF(a s)) x g 


i=0 





x=0 

pa 


d xe~ xVTsF{as) q d { A= 



— G(y = as) . 
7 T a 


Note that in the discrete sum over x in the first line of Eq. (4.6) only the terms such 
that x + a is even contribute - while the terms such that x + a is odd are just zero. Hence 
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this discrete sum is approximated by 1/2 times the integral over x , as indicated in the 


second line of Eq. (4.6). The function G(y) in the third line of Eq. (4.6) is given by 

G(y) = 1 - y/FyF(y) exp [yF 2 (y)]evic [y/yF(y)\ . (4.7) 

Finally, one obtains the asymptotic behavior of the generating function of Q%(n)( o) in 


Eq. (4.2), as z —> 1 




Q 


n =0 


y/l-z ' CQ V vr J 0 yjy 

such that Q%{n )( 0 ) ~ 2 cfk / yjn r n as n —> oo, for n even (while Q B (n)^ — 0 if n is odd). 
Using that G d (0,0,n) = 2 _n ( n " 2 ) ~ ^2/7rn~ 1 ^ 2 , for n » 1, one obtains from Eq. (4.1) 
together with Eq. (|4~8) 


^-•Gfe), 


(4.8) 


2 7°° d v 

lim Q%(n) = Q%( oo) = V2c% = -= / — e~ y G(y) = 0.6543037... , (4.9) 

rwoo v V 71 Jo Vv 

which is different from the corresponding value for free random walks, Q d (oo) = 


0.626508... [see Eq. (1.17)]. 


4-1.2. Exponential jump distribution We now turn to the computation of Q%(n ) in 


Eq. (1.7) in the case of an exponential jump distribution. In this case, our starting 


point is the joint distribution of the ages of the records given in Eq. (3.22). As before 
in the case of discrete random walks (|4.1[), one has 


W = X] 5Z m ’ n ) ’ 

m =1 a =1 £ 1 = 1 


(4.10) 


from which, using Eqs. (3.22) and (3.30), one obtains [similarly to the discrete case in 

Qb( u )( o) 


Eq. (4.2) 


Q e B (n ) = 


OO 


n =0 


G e (0,0, n) 

OO OO 

<%(")(0)2" = 27*° 27lV(s,a)] m “V(m.o) , 


(4.11) 


a=1 


m= 1 


where f e (i) is given in Eq. ( 3.27[ ). 

To obtain the large n behavior of Q e B {n)^p we need to analyze its generating 


function in Eq. (4.11) in the limit z = e s —> 1, i.e. s —> 0. In this limit, the sum over 


a is dominated by large values of a ~ 1/s. In this limit, h e (z : a) takes the scaling form, 
for s —)• 0, as fixed 

h e (z, a) ~ (l — y/sF(a s)) , (4-12) 


where the scaling function F(y) is given in Eq. (4.5). Furthermore, in the limit 


z = e s —» 1 (i.e., s —> 0), the sum over m in (4.11) is dominated by large values 
of m. In the limit, m —> oo, a —>■ oo, keeping m/y/a fixed, we obtain, from the GF of 


g e (m, a) given in Eq. (3.32), that it takes the form 

e / , 1 

q (m, a) 


2 b^J~T\ a 3 / 2 


m c -m 2 /(4a) 


(4.13) 
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Therefore, by injecting these asymptotic behaviors given in (4.12) and in (4.13) into Eq. 
( TO , one obtains 




sa 




171 ^ —m 2 /(4a) 


(4.14) 


2 b^/n „ 

n =0 a=l m =1 

The sum over m can then be approximated by an integral, in the limit s —» 0 as follows 


y m (l — \fsF(a s))e 


m 1 —m 2 /(4a) 


m= 1 


E 

m=l 


m g-(m-l). sfsF(as) e -m 2 /(Aa) 


d m me- <m - ,, ^ F( ‘“>- m2/ < 4a> = 2aG(as) , 


(4.15) 


where the function G(y) is defined in Eq. (4.7). Hence, finally, approximating the 


remaining sum over a by an integral (in the limit s —» 0) and performing the change of 
variable y = sa one obtains 




-Q 


3-2/ 


n=0 


\/l — -2 


i C Q — 


T—= / dy—G(y) , 

Jo Vv 


(4.16) 


which implies that, for large n, Q e B ( n )( o) ~ c q/V™- Hence, using that G e (0,0,n) ~ 
l/(26y / 7m), for large n, one obtains from the first line of Eq. (4.11) that 

hm Q%(n) = Q e B ( oo) = Q%( oo) , (4.17) 

n—¥ oo 

which is thus the same constant as the one appearing in the discrete case given in 


Eq. (4.9). 


4-2. Age of the longest lasting record 


4-2.1. Discrete random walk. The average value of &xB (n)) defined in Eq. (1.8) can 
be computed from its cumulative distribution as follows 

OO 

(Cx,bM} = Eh - F d (e,n)) , F d (t,n) = Pr[C„, B (n) < f) , (4.18) 


e=o 


where F d (£,n ) = Pr[f ^, av R (n) < 1} is simply obtained by summing up the joint PDF 


m 


Eq. (|3.7|) over £i, ■ ■ ■ , £ m _\ and a from 1 to i - as all of them have to be smaller than i 

ollows: 

i 

y P d {£- ,£m-i,a,m,n) .(4.19) 


- and finally over the number of records m as follows: 

oo i l £ 

F\l,n) = Pr [Cx „M < t] = E E E ■ 


m=la=l£i=l 


Therefore, one has 


(Cax,sW) = 


(t 


max,B 


n 


))( 0 ) 


G d (0, 0, n) 
where the numerator is given by 

OO t 


(4.20) 


<4 


max, B 


»)(o) = G d (0,0,n) - yyy ■■■ y P d (£ 1 ,£ 2 ,--- ,£m-i,a,m,n) {0) 


m=l a =1 £\=1 


(4.21) 
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CL-\-\ 


bM}(o>v = E g“(o. o.*) - E E EW)] G l(™ - ro.“) 


n=0 


£=0 


a=l m=l 


(4.22) 

where the function h d (z,£) is defined in Eq. (4.3). In the limit z = e _s —> 1, i.e. s —> 0, 


one can replace G>(m — 1,0, a) by its scaling form given in Eq. (4.4) and h d (z ,£) by 
its asymptotic behavior given in Eq. (4.5), in the scaling regime where a and £ are 


large, keeping the products sa and s£ fixed. The sum over m in Eq. (4.22) can then be 


analyzed along the same lines as done before in Eq. (4.6) to yield: 

CL~\~ 1 

^2[h d (z,£)] m 1 G>(m — 1, 0, a) 


m= 1 


~ \I~~ [l — saF(s£) exp [saF 2 (sf)]erfc [-\/saF(s^)]] , (4.23) 

where the function F(y) is defined in Eq. ( |4.5[ ). Inserting this asymptotic behavior 
(4.23) in Eq. ( |4.22[ ), one finds 


£<£«»>(o)*" ~ 


, as z —> 1 , 


(4.24) 


where the amplitude c^ ax is given by [remembering that only the even values of n 

(4.25) 


contribute to the left hand side Eq. (4.24)] 


r c 

= V2 

Jo 


dt 


2 - I(t) 


where the function I (t) is given by 

1 /' t d3/___ 

'o Vv 


m = ~,= f (l - \J~FyF{t)e yp2 ^erfc (y/yF(t)j) 


F(t)e t+tF ^erfc[-\ /tF(t)] — e 


(4.26) 


(4.27) 


1 — F 2 (t) 

The behavior in Eq. (4.25) implies that /maxsl 71 )) ~ (4 c max/ a/tt) y/n, for large n (even). 


Using finally that G d (0, 0, n) rv_/ a/2/ (nn) (again, for n even) one finally obtains 


lim 

n—> oo 

= 4 


(Cax,/ W )) 

n 

coo / 


_ \d 

/ 'max,i3 


= 2/2 < 


dt ( - — J(t) 


= 0.6380640 


(4.28) 

(4.29) 

(4.30) 


which is strictly smaller than the constant Qf 1 ( 00 ) given in Eq. (4.9) 


4-2.2. Exponential jump distribution. The average value of xg (n)) defined in 


Eq. (1.8) can be computed from its cumulative distribution as follows 

OO 

<C* B (n)> = E* 1 - f ’(<») - F ‘(en) = Pr[C«, B (») < <\ . (4.31) 


e =o 
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where F e (l,n) = Pr[^ nax B {n) < £\ is simply obtained by summing up the joint PDF in 
Eq. (3.22) over l\, - ■ ■ ,l m -i and a from 1 to l and finally over the number of records 
m as follows: 

oo l £ £ 

F‘(l,n) = Pr[C», B («) < «] = EE ZD' E **(*1. V" ,C-i,o,m,n) . (4.32) 

m =1 a=l £i=l Cm- 1=1 

Hence (l e mSb xB (n)) can be written as 


(P 


e 

//e / \\ '~max,jj\' ~n W 

Ua x,M P G e (0,0, n) ’ 


4>(o) 


(4.33) 




max,!? 


»)( 0 ) = 


£=0 


oo £ £ £ 

G e (0,0,n) - pe (^ir-- ,£m-i,a,m,n) {0) 

m=1 a =1 q = l 4rt-i=l 


Again, to extract the large n behavior of b ( u ))( o ), it is convenient to analyze the 
generating function of (4ax,.B ( n )) (o) > which can be easily done by using the previous 
analysis. Indeed, one has 


G e (0, 0, e -s ) = ^G e (0,0,n)e 


n =1 


~ sn ~ \-= , ass^O 
2b ^ ’ 


(4.34) 


On the other hand, from the analysis performed above in Eqs. (4.12, 4.13, 4.15) one 
has 

oo oo £ i 

^ e ~ Sn ^2^2" ' pe (^i>'" ,tm-i,a,m,n) (4.35) 

n =1 m=l £i=l £ m - i=l 

[l — y/nas F(s l) exp [asF 2 (s£)]erfc [y/saF(s£)]] , (4.36) 


in the limit s —» 0, a —> cxo, £ —> oo, keeping sa and si fixed and where the function 
F(y) is defined in Eq. (4.5). Therefore, using Eqs. (4.34) and (4.35), we obtain 


r e 

E<W")>(o)z”~ b{l ”“ 3/2 


, as z —> 1 , 


where the amplitude c^ ax is given by 


max 

oo 


d*( 2 " ) ’ 


(4.37) 


(4.38) 


where the function I(t) is given in Eq. (4.26). From Eq. (4.37), one deduces that 
(Cax,BW>(o) ~ (2 c e m J[b for large n. Using that G e (0,0,n) ~ 1/(2 by/wn), for 

n 1, we obtain finally 


,. (4a X ,bH) \ e A e 

llm - = 4ax,B = 4 C max 

n—>oo n 


= 4 /o dt l2- /(t) 

= 0.6380640 


(4.39) 

(4.40) 

(4.41) 


which is exactly the same constant that appears in the discrete case in Eq. (4.30). 
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We conclude this section by the study of the full probability distribution of 
t“ a x b( n )■ We perforin here the analysis on the exponential case, but it can also be 
done on the discrete case, leading to the same result in the large n limit. For the 
exponential case, one obtains the GF of F e (£, n)( 0 ) = Pr[l ^ a Y R (n) < /]( 0 ) from Eqs. 
(4.33) and ( |3.22[ ) as 

i 


OO 

£ 

n =0 


F"(£,n) (0)2 “ = £>“£;[ 


h e (z,£) 


I m —1 


a =1 


m =1 


q (ra, a) 


(4.42) 


Again, to investigate the large n behavior of F e (£,n)( o), we need to analyze this GF 
in the limit 0 = e -s —y 1, i.e., s —> 0. From the analysis performed before, using in 


particular the asymptotic behaviors in Eqs. (4.12) and (4.13) one obtains 


OO 


ti b ^ s 


(4.43) 


where the function /(f) is given in Eq. (4.26). From Eq. (4.43), we obtain that, for large 
n, F e (£,n) is a function of the ratio £/n, as expected from our previous computation 
showing that R (n)) ~ R n (4.41). This expresses the fact that, for large n, 
the random variable lZ(n) = 4axs( n )/ n becomes independent of n, lZ(n) —> 1Z when 
n —> oo, as it is also the case for free RWs [29], The cumulative distribution of this 
random variable R, F R (r) is given by 


F’(t,n)=F^ - = r 


n 


^(r) = 2 


/ (u r) 


u 


(4.44) 


where T is a Bromwich contour in the complex plane and r G [0,1]. From this expression 


(4.44), it is possible to study various features of the PDF of R, /^(r) = dF^(r)/dr. 
In particular, one finds that /^(r) is non-analytic at r = 1/2 and one can compute it 
explicitly on the interval [1/2,1] where it takes a simple form 


fn( r ) = 1 + 


4 r 3/2 


, - < r < 1 
’ 2 ~ ~ 


(4.45) 


On the other hand, for small r, it has an essential singularity where it behaves as 


(4.46) 


In fn( r ) ~ ^ , as r 0 . 

The study of the corresponding PDF for free RW was recently carried out in Ref. 

(see also Ref. [30]). In this case, the PDF also exhibits a non-analyticity at r = 1/2 
but in this case, it has a square root singularity when r —> 1, while for a bridge, 
there is no singularity [see Eq. (4.45)]. On the other hand, the PDF for free RW also 


exhibits an essential singularity when r —>• 0 as in Eq. (4.46) but without the logarithmic 
correction, which is thus a special feature of the RW bridge. In the next section, we 
present a numerical evaluation of this PDF /^(r) (see the left panel of Fig. [TT|. 


5. Numerical results 


In this section, we present results obtained from numerical simulations of random walk 
bridges, which we compare to our exact analytical computations. We consider discrete 
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(i a = d) as well as continuous (a = c) jump distributions, for g. = 2, with a special 
attention on the exponential jump distribution (a = e). 


5.1. Discrete random walk bridges 


To simulate discrete RW bridges, we generate simple free random walks of n steps 


according to Eqs. (1.1) and (1.2) and retain only the walks such that x{n ) = x(0) = 0. 
The statistics is then performed over this restricted ensemble of RW bridges. In Fig. [3j 
we show a plot of the average number of records (i?^(n)), computed numerically, as 
a function of yjn , which shows a perfect agreement with our exact results in (2.10). 


In Fig. [3] we have also plotted the asymptotic estimate of (R%(n)) beyond the leading 
0(y/n) term. Indeed, from the exact formula in Eq. (2.10), one obtains: 


(R 


B\ n )) 


2 3/2 


n +2 + ’ 


(5.1) 


which, as can be seen on Fig. [3J is a very good estimate of the exact result for y/n > 5. 



Figure 3. Average number of records (R%(n)} as a function of y/n for discrete RW 
bridges. The full blue circles are the results of numerical simulations. The green empty 


circles correspond to the exact results given in Eq. (2.101, while the full red line is a 


plot of its asymptotic form given in Eq. (5.11. For comparison, we have also plotted, 
with empty blue circles the exact results for ( R d (n )) for the free RW, as obtained in 
[5], which displays a faster growth (R d (n)} ~ (2/ y/n) y/n for n^> 1. 


We have also computed numerically the distribution of R^(n), 


Pg(m,n) = 


Pr (Rg(ri) = m) both exactly, for n — 40 and 80 from the analytical formula in Eq. (3.13) 
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and numerically for n = 40. As shown in Fig. [4| the numerical results agree perfectly 
with our exact analytical results. Besides, for large n and large m, keeping m/y/n = X 



Figure 4. Plot of \fn T 3 / (m, n) as a function of (m — l/2)/-y/n for discrete RW - 
see Eq. (5.2) and the discussion below in the text. The filled circles correspond to 
the exact values obtained from Eq. (3.13) for different values of n = 40 (blue) and 
n = 80 (green). The open circles correspond to the results of numerical simulations 
for n = 40, which coincide exactly with the exact results. Finally, the solid line is the 
exact scaling function ipg(X) given in Eq. (1.19). 


fixed, one expects the scaling form given in Eq. (1.20). Our simulations for finite values 
of n exhibit small finite size corrections to this limiting scaling form. As suggested 
in Fig. [4| these finite n actually correspond to a simple shift of the scaling variable 
m/y/n —* (m — 1/2)/a /n, i.e., 


\ n 


m - 1/2 \ 

n J 


(5.2) 


where the scaling function (p] 3 (X ) is given in Eq. (1.20). Although we have not 
performed a detailed analytical study of the finite n corrections to the limiting PDF 
in Eq. (1.20), this form in Eq. (5.2), and in particular this shift of 1/2, is consistent 
with the asymptotic expansion of the average number (Rg(n)), beyond the leading 
<D(y/n) term, given in Eq. (5.2). 

Then, we have computed numerically the probability of record breaking Q%(n), 
both from numerical simulations and from our exact formula in Eq. (4.2). As shown in 
the left panel of Fig. [5j the agreement between both estimates is very good. Besides, 
we also see on that figure that the convergence of Q%(n) to its asymptotic value oo) 
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is quite slow. The right panel of Fig. [5] indicates that the speed of convergence is 
proportional to 1/y/n, i.e., Q%(n) — Q B (oo) oc l/y/n. Furthermore, the estimate of 
Qb( oo) obtained in this way, by extrapolation, is in good agreement with our exact 
results in Eq. (4.9), see the right panel of Fig. [5j 



l/y/n 


Figure 5. Left: Probability of record breaking Q%(n) for the discrete random 
walk bridge. The blue open circles correspond to our numerical simulations while 
the yellow full circles are the exact values of Q B (n), extracted from the generating 
function in Eq. (4.2). Right: Slow convergence of Q B (n) to its asymptotic value 


Q'g(oo) = 0.6543037... in Eq. (4.9). The full line is a guide to the eyes and corresponds 
to Q%( oo) + const. /y/n. 


Finally, in Fig. |6j we show numerical results for (^ a X)B (n))/n, obtained 
from numerical simulations and compare them to the exact results which can be 


straightforwardly obtained from the generating function in Eq. (4.22). Here also the 


agreement between numerics and theory is very good. We also notice that, as in the 
case of the probability of record breaking (see Fig. |5j), the convergence of (^ aX)S (n))/n 
to the asymptotic value A^ ax B 


see Eq. (1.25)] is actually quite slow. 


5.2. Continuous jump distribution with /i = 2. 

To generate a RW bridge with continuous jump distributions and p, — 2 (i.e., with a 
finite variance cr), we first generate a free random walk of n steps, starting from x(0) = 0 


and evolving according to Eq. (|l.l|). We subsequently use the following construction 

(5.3) 


k 


x B (k) = x(k) - x(n) 


n 


to generate a RW bridge {x B (k)}, for 0 < k < n. In particular, Eq. (5.3) implies 


obviously x B (0) = x B (n) = 0. This simple construction (5.3) holds, for any finite n, for 


Gaussian jump distributions. For other jump distributions, including in particular the 
exponential jump distribution (a = e), this construction is only approximate, becoming 
exact in the limit of large n, where x B (k)/y/n, for k ~ 0(n), converges to the Brownian 
bridge. In Fig. [Tj we show a plot of the average number of records ( R B (n )) both 
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Figure 6. Plot of B (n))/n as a function of n for the discrete random walk bridge. 
The open circles correspond to the exact analytical value computed from the generating 


function in Eq. (4.22) while the full circles correspond to the estimates obtained from 
direct numerical simulations. The solid line corresponds to the asymptotic value 
max.it = 0.6380640..., see Eq. (1.25). 


for an exponential distribution, a = e, and for a Gaussian distribution, both of them 


belonging to the class of continuous jump distributions a = c as in Eq. (1.3) with /i = 2. 
For large n, one expects that, in both cases, ( R%(n )) rv_/ (R c bH) ~ a b(r = 2 )V^ 


with A c B (n = 2) = yTr/2 [see Eqs. (1.21) and (1.22)], in good agreement with our 
numerical simulations (see Fig. [7J). In fact, our numerical simulations suggest that 
the first correction to this asymptotic behavior is the same for the exponential and the 


Gaussian cases and can be computed from our exact formula for ( R e B (n )) in Eq. (2.17), 
yielding 

1 + 0(1). (5.4) 


Note that the first 0(1) correction, namely 1/2, is the same as in the discrete case (5.1). 
For smaller values of n, the data for the exponential case, a = e, in Fig. [7] show a slight 
discrepancy between the numerical results and our exact formula for ( R e B (n )) given in 
Eq. (2.17). This is due to the fact that our construction of the RW walk bridge ( |5.3 ) is 
in this case only approximate for small values of n. Note that to simulate exact bridges 
in this exponential case, one could use the Monte-Carlo method used in Ref. ra. which 
however requires more numerical efforts and goes beyond the scope of the present work. 

In both cases, exponential and Gaussian jump distributions, we have computed the 
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Figure 7. Average number of records {R%(n)} obtained by simulations of a RW 
with an exponential jump distribution, a = e, (blue circles) and for a Gaussian jump 
distribution, a = c, (green circles), compared to the exact asymptotic value (full 
black line), given in Eq. (5.4). The red circles correspond to the result for the 
exponential case in Eq. (2.17). The slight discrepancy between the exact and the 
numerical results in the exponential case for small values of n is due to the fact that 
our numerical procedure to generate a bridge ( |5.3| ) is only approximate for finite n. For 
comparison, we have also plotted (with cyan circles) the exact results for a free RW 
with a continuous jump distribution obtained in jS], which displays a faster growth 
( R c (n )) ~ 2/y/ny/n for n^> 1. 


distribution of the number of records P B (m,n) and P B (m,n). The results are shown 
in Fig. [8j In the left panel, we show the results for an exponential jump distribution, 
a = e, which show a very good agreement between our numerical estimates and the exact 
formula for P B (m,n) given in Eq. (3.34). These data for finite n show some deviation 
from the asymptotic result in Eq. (1.23) valid for n » 1. As suggested in Fig. |4j and 
similarly to the discrete case (5.2), these finite n corrections actually correspond to a 
simple shift of the scaling variable m/y/n -A (m — 1/2 )/y/n, i.e., 


P e B {m,n) - 


n 


m 


FA 

n ) 


(5.5) 


where the scaling function ip e B (X) is correctly predicted by the expression given in 
Eq. (1.23). Note that this shift of 1/2 is fully consistent with the asymptotic expansion 
of (R B (n)) beyond the leading order in Eq. (5.4). In the right panel of Fig. [8j we 
show numerical data for P B (m,n ) for a Gaussian jump distribution. These data are 
particularly instructive, as we do not have any analytical result for this quantity for 
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Figure 8. Left: Plot of \fnP B (m, n) as a function of (to —l/2)/y / n for an exponential 
jump distribution. The filled circles correspond to simulations for n = 40 (blue circles) 
and n = 80 (green circles). The empty circles correspond to the exact analytical result 
given in Eq. ( |3.34 1 for n = 40. The solid line corresponds to the exact scaling function 
(p e B {X) given in Eq. (1.231. Right: Plot of y/nP B {in., n) as a function of (to — l/2)/- v /n 
for a Gaussian jump distribution for n = 160 and n = 320. The solid line corresponds 
to the exact scaling function gf B (X) given in Eq. (1.23) 


Gaussian jumps. Our data actually suggest that the distribution of the number of 
records in this case is also described, for large n, by the same scaling form as for the 
exponential jump distribution in Eq. (5.5), with the same scaling function (p B (X). These 
data support our conjecture that this result in Eq. (1.23), which we can explicitly show 
only for the exponential case, should actually hold for any continuous jump distribution 
as in Eq. (1.3) with p = 2. 

We have also computed numerically the probability of record breaking Q e B (n ) for 
an exponential jump distribution, which we have plotted in Fig. [9j Here again, we see a 
very good agreement between our exact results [extracted from the generating function 
in Eq. (4.11)] and the numerical estimates. These data also corroborate the asymptotic 
behavior Q%(n ) —>• Q B (oo) = 0.654304... when n —* oo, where Q B ( oo) = Q B ( oo) is 
given in Eq. (4.9) and coincides with the asymptotic value Qf j(oo) found for the discrete 
RW. Note that our numerical results indicate that the convergence to the asymptotic 
value is oc 1/n, which is thus faster than for the discrete case, where the convergence is 
1 /a Jn (see the right panel of Fig. [5|. In Fig. [9j we have also plotted the estimates 


oc 


of Q c B {n ) for random walk bridges, with Gaussian jump distribution, obtained from 
numerical simulations, using Eq. (5.3). The data for exponential and Gaussian jumps 
are, for n sufficiently large, indistinguishable, which reinforce our claim that all our 
asymptotic results for the exponential case should also hold for any continuous jump 
distributions with // = 2, in the limit n —> oo. 

we show a plot of B (n))/n, for an exponential jump 


Similarly, in Fig. 


10 


distribution and of B (n))/ri, for a Gaussian jump distribution. In the exponential 
case, we have plotted both our exact results, which can be derived from Eq. (4.33), and 
the estimates obtained by simulating the random walk bridge, using Eq. (5.3). We see 
that the agreement is very good. Interestingly, we also see that the data for exponential 
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Figure 9. Probability of record breaking for an exponential jump distribution, Q e B (n ), 
and for a Gaussian jump distribution, Q B (n). In the exponential case, we show 
data corresponding to the exact result (open circle), obtained from the generating 
function in Eq. (4.11) and data obtained from direct numerical simulations (full circles) 
using Eq. (5.3). The data for a Gaussian jump distribution correspond to numerical 
simulations using Eq. (5.3). The convergence towards the exact asymptotic value 
Q%( oo) = 0.654304... is fast (ex 1/n). 


and Gaussian jump distributions are almost indistinguishable for n sufficiently large. 

Finally, we present the numerical evaluation of the distribution //j/n) = dF^(r)/dr 
°f Caxs( n )/ n i 11 large n limit [see Eq. (4.44)]. The numerical estimates of /^(r) 
can be obtained by evaluating F e (£, n) from the generating function that was explicitly 
computed in Eq. (4.42), which we expand close to z = 0. The same procedure was used 
by us in Ref. [22j to plot the corresponding PDF for the free RW. The plot of /^(r), for 
RW bridges with an exponential jump distribution, obtained this way is shown in the 
left panel of Fig. 11 The non-analyticity at r = 1/2 is clearly visible and in addition, 
we observe a perfect agreement between the numerical evaluation of /^(r) and the exact 
result for 1/2 < r < 1 [see Eq. (4.45)]. Finally, in the right panel of Fig. 11, we show, 
for completeness, a plot of the PDF of V = 1/TZ, whose PDF turns out to be simpler to 
study (see Ref. [2.9] as well as Appendix C). 


We conclude this section by noticing that a natural way of simulating numerically 
random walk bridges would be to write an effective local equation of motion for xs(k) 
which would take into account the conditioning that the walker has to return back to 
the origin after n time steps. For Brownian motion, which is continuous both in space 
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Figure 10. Plot of B (n))/n, for an exponential jump distribution, and 

(4ax b ( n))/n, for a Gaussian jump distribution, as a function of n. For the exponential 


case, the exact results are obtained from Eq. (4.33), while in both cases, the numerical 
estimates are obtained by simulating the random walk bridges from Eq. (5.3), with 
the appropriate jump distribution. 




Figure 11. Left: Probability density of the ratio TZ(n) = ^ max (n)/n for an exponential 
jump distribution, for n = 30 and n = 40. The points were obtained from the exact 
expression for the generating function of the cumulative distribution F e (£,n) (4.42). 
The red curve is the analytical prediction for the limiting PDF /Jj(r) (4.45) valid for 
1/2 < r < 1. Right: Probability density of V(n) = n/£ max (n) for an exponential 
jump distribution of steps, for n = 30 and n = 40, obtained from the left panel. 


and time, this can indeed be done by writing an effective Langevin equation for the 
Brownian bridge [33 . 34f[ . Extending this approach to discrete time random walks - 
both with discrete and continuous jumps - is a very interesting open problem. 
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To conclude, we have studied the record statistics of RW bridges, for different types 
of symmetric jump distributions. Our results show that the record statistics of such 
constrained random walks, which are constrained to start and end at the origin after 
n steps, are quantitatively different from their counterpart for free RWs, which can 
end up anywhere on the real axis. We first showed that the statistics of the number 
of records R%{n) is not only different for discrete (a = d) and continuous (a = c) 
distributions but, even for continuous jump distributions, also depends on the details 
of this distribution. We obtained exact results for the average number of records 
(«SW) rv_/ A%y/n and showed in particular that A C B = A B (/u ) depends continuously 
on the Levy index /j characterizing the RW. This is quite different from the free RW 
case where R c (n ) ~ A c y/n where A c = 2 /^/tt, independently of /i. Furthermore, we have 
computed exactly the full record statistics for two different types of jump distributions: 
the discrete RW (a = d) and the exponential distribution (a = e). We emphasize that 
these calculations are technically much harder than for free RWs as the records statistics 
of RW bridges requires to keep track not only of the ages of the records but also of the 
increments between successive records [see Fig. [I] and Eq. (3.21)]. For these two jump 
distributions, we have also computed the probability of record breaking Q B (ri) and the 
average age of the longest lasting record (f" ax B {n)) and shown that they give rise to 
two new non-trivial constants which we have computed explicitly (see Table 1). 

Although the records statistics for a continuous jump distribution depends, for finite 
n, on the details of this distribution, one expects that, for large n, it only depends on 


the Levy index /i in Eq. (1.3). Although we can not prove this statement, our numerical 


data indicate that this should indeed be the case, at least for p = 2. This implies that 
our asymptotic results (see Table [Tj) obtained in the exponential case, which we could 
solve exactly, should describe the large n limit of any continuous jump distribution with 
fi = 2. The generalization of our results to arbitrary value of the Levy index 0 < p < 2 
remains a challenging open question. 
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Appendix A. Useful formulas for discrete random walks 

Here we consider a discrete RW, starting from x(0) = 0 and evolving via 

x(k) = x(k — 1) + r](k) , (A.l) 

where the jump variables rj(ky s are i.i.d. random variables distributed according to 


p d (rf) = \S(rj +1) + lS(rj — 1). Let us denote by W(x, n ) the number of lattice RW (A.l) 
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starting at x(0) = 0 and ending in x after n time steps. One has 

(n+x) , if n + x is even 
V 2 7 

0 if n + x is odd . 

From W(x,n), one obtains immediately the discrete propagator G d (x,x 0 ,n ) as 


W(x,n ) = 


35 


(A.2) 


-.a. \ W{x-x Q ,n)\ h( n+x-xp ) , if n + x- x 0 is even 

G d (x,x 0 ,n)= 1 • ,, ( A - 3 ) 

z I (J it n + x — x 0 is odd . 

The generating function G d (x, xq,7i) of G d (x, xo,n) can be computed from this explicit 
expression (A.3) and it yields: 


G d (x,x 0 ,n ) = ^G d O,a;o, 

n=0 


n)z'‘ = 


E 

fc=r^i 


z ^ 2 k-(x-xo) f2k — (x — x 0 ) 
k 


vT 


i - 


X—XQ 


, (A.4) 


where [u] is the smallest integer not less than u and where the last equality can be 
obtained using, for instance, Mathematica. 

Furthermore, from the expression of G d (x,xo,n) (A.2), one can also compute the 
constrained propagator G>(x,0,?z) using the method of images. One has indeed: 

G>(x, 0, n) = G d fx, 0, n ) — G d (x, —2, n) (A.5) 

A — (n±^ +1 ) j , if n + x is even 

0 , if n + x is odd 

The generating function G>(x, 0, z ) of G%{x, 0, n) can then be obtained from this explicit 
expression (A.5) as 


(A.6) 


G>(x,0,z) = G>(x,0,n)z r ‘ 


n =0 
oo 

E 

fc=rn 


2 (1 - Vl - 


z \ 2 k—x (2k — x 
k 

x+l 


E 

fc=r^i 


z ^2k-(x+2) {2k- (x + 2) 

k 


(A.7) 


Note that this result (A.7) which we obtained here by counting paths can also be 
obtained using a backward Fokker-Planck equation. 

Finally, one can also compute the propagator G d (x, 0, n) associated to a RW which 
is constrained to stay strictly above 0. To do so, one can simply analyze the first 
step. This first step is necessarily a step +1, such that rj(l) = +1, which happens 
with probability 1/2. After this first step, the remaining propagator is just given by 
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— 1, 0, n — 1). Therefore, one has 
G>(x, 0, 0) = < 5.^0 

G > (x, 0, n) = -G>(x — 1, 0, n — 1) , n > 1 . 
One thus obtains the generating function G > (x , 0, z ) as 


G>(x,0,«) = J^G>(a;,0,n)z n = 

71=0 


VT 


(A. 8 ) 

(A.9) 

(A. 10 ) 


Appendix B. Useful formulas for RWs with exponential jump distribution 

Here we consider a RW, starting from x(O) = 0 and evolving via 

x(k) = x(k — 1) + rj(k) , (B.l) 

where the jump variables 77 (h)’s are i.i.d. random variables distributed according to 
P e (v) = 2 ^e _ l r? l //b . A remarkable feature of the exponential jump distribution is that the 
constrained Green’s function G^x, x 0 , n) can be computed exactly [IHTj . 

To begin with, it is useful to compute the Fourier transform of the jump distribution 

1 


from which we easily obtain the free propagator 

(Tf r r. n I = / _ - _ 

' ’ [l + (6?) 2 ] 

Finally, the GF of G e (x, xo, n ) is given by 

OO 

G e (x, x 0 , z) = G e (x, x 0 , n) 


z n = z 


n =1 


d q 
2 tt ' 


iq(x-x 0 ) 


1 — z + (bq 2 ) 


(B-2) 


(B-3) 


(B-4) 


where the Dirac delta function S(x) comes from the term n = 0, while the second one 
comes from the sum from n — 1 to oo of the geometric series. Finally, the integral over 


q in (B.4) can be explicitly evaluated to yield the result given in Eq. (|2.12|) 

OO 

z 


G e (x, 0, z) = ^VG e (x,0, 


z = 


| a: | n - 


(B.5) 


2by/l — z 

n= 1 

We now come to the constrained propagator Gf(x, 0,n), which is harder to 
compute. The easiest way to compute it is to use the so-called Hopf-Ivanov formula 
[ 35] (see also [36] for a detailed derivation of this result) which gives the following 
expression 


E-~7 

71=0 

where 0 (A, z) is given by 


(f(X,z) = exp ( -- 


d q 


; - a " = 0(A,^), 

(B. 6 ) 

In (1 — zp e (q))\ 

X 2 + q 2 )' 

(B.7) 
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where p e (q) is the Fourier transform of the jump distribution given in (B.2). The function 
( p(X,z ) is thus given by 


0 (A, z) = exp 


A 


7 r 


dg 


1 ( 1 -z+(bq ) 2 

111 ^ 1 + (&g) 2 

A 2 + q 2 


(B.8) 


The integral over q in Eq. (B.8) can be explicitly evaluated, using the formula 4.295-7 
of Ref. [37] 


00 ln(a 2 + 6 2 a; 2 ) n fad+ bc\ 
da: —--; In { --- | , tor a,b,c,d > 0 , 


cd \ d 
1 +A b 


(B.9) 


j 0 c 2 + d 2 x 2 

to yield finally the identity 

OO /»00 

dxG e > (x,0,n)e~ Xx = 

n =0 

It is then easy to inverse the Laplace transform with respect to x to obtain Gt,(x,Q,z) 


y/l — z + A b 


= 1 + 


1 - v^ 


\Jl — z + A b 


(B.10) 


given in Eq. (2.13) 


Gl(x, 0 ,z) = J2 z n Gl(x, 0, z) = S(x) + 


1-VT 


Z _ M 

—e b 


(B.ll) 


n=0 


Note finally that, by differentiating the relation in Eq. (B.10) (m — 2) times with respect 


to A and setting A = 1/6, one obtains straightforwardly the relation given in Eq. (3.32) 


Appendix C. Analysis of the distribution of l ]/ a Y R (n)/n 

In this appendix, we derive the behavior for the PDF f%fr) first for r G [1/2,1] and 


then for r —» 0 given in the text in Eqs. (4.45) and (4.46) respectively. 


Appendix C.l. Density f^{r) in the interval [1/2,1] 


We start from Eq. (4.43[) given in the text 


Y_ e - n F‘(e,n) lm = , 


n>0 


(C.l) 


valid in the limit s —)• 0, £ —)• oo, keeping the product x = s£ fixed and where the 


function /(f) is given in Eq. (4.26). In the limit s —> 0, the discrete sum in Eq. (C.l) 
can be replaced by an integral and F e (/, n)( 0 ) can be written as 




(C.2) 


F e {£,n) { o) = Prob(Ca*,fl(w) < ^)(o) = Prob 


where T is a Bromwich contour. Note that by definition, we have 

^ >?),„,- < a3 > 
In the limit of large n this probability converges to the complementary distribution 
function of the random variable V — 1/7Z = lim^oo n// max (n), 

(C.4) 


= Prob (V > v) m 
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with v = n/i and x = s£ as above. Finally, dividing by G e (0,0,n) ~ l/(2by/im), 

_ _ f r\x 1 

KM = J (°-5) 

We want to analyze the behavior of F^(r) for r close to 1, or equivalently the small v 
behavior of of K(v), which is dominated by the large x behavior of the integrand in 


the Bromwich integral in Eq. (C.5). Tims we need to expand I(x) given by (4.26) in 


the small parameter e(x) = e x . Formally 
P Hr 1 f nr 1 

-er-pW c) « / — (l 0 (x) + X,(x) + X 2 (x) e^ + ...), (C.6) 

W X ./ 2l7T Jx 


2m \fx J 2i7r 

where the functions Tf.{x) admits an expansion in powers (both positive and negative) 
of y/x. Therefore the first term oc e~ x contributes to the function for v > 1, the second 
term, oc e~ 2x , for v > 2, and so on. At first order in e(x) = e~ x we have 




— {Ax + l)exicy/x J + 0{e 2x ) , 


which, divided by y/x, can be Laplace inverted to give for F v (v) 


W = ) + i-A (1 < u < 2) , 


and therefore 

f» = -4-f^w = 4 + tA 

av v z Ay/v 

Finally, using f%{r) = /y(l/r)/r 2 , one obtains 


(1 < v < 2) . 


fn( r ) = 1 + 


4 r 3/2 


- < r < 1 
2 




(C.7) 

(C.8) 

(C.9) 

(C.10) 


as announced in the text in Eq. (4.45) 


Appendix C.2. Small r behavior of f^{r) 

We now study the small r behavior of /^(r), or equivalently the large v behavior of 
fv(v) (recalling that V = l/7£). One can show that the function I{x)/yfx is an entire in 


the complex x plane. Hence the integral (C.5) is expected to be dominated by a saddle 
point in x. Let us define 

I{x) 


4>{x) = In 


x 


(C.ll) 


The saddle-point equation is v + 4>'{x) = 0, thus (p'{x) is large and negative, hence x is 
large and negative too. Setting x = —y 2 we obtain the estimate the large negative x 


behavior of F(x) in Eq. (4.5) as 


which yields 


F(x = - y 2 ) 


<t>{x = -y 2 ) 


2yfi n/ 3 ’ 






Airy 4 Attx 2 ’ 


(C.12) 


(C.13) 
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coming from the dominant term ein I(x). The saddle-point equation now gives 
v ~ e~ 2x , (C.14) 

hence, for v large, using estimates at exponential order, we find 


KM = ~ 


(C-15) 


i.e., a super-exponential decay for this quantity. As a consequence one can also predict 
that the essential singularity at the origin of /^(r) = /y(l/r)/r 2 is given by 


/i?M 


3 lnr/(2r) 


(C.16) 


as given in the text in Eq. (4.45). 
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